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2 Temporal behavior of quantum mechanical systems 



The temporal behavior of quantum mechanical systems is reviewed. We mainly focus 
our attention on the time development of the so-called "survival" probability of those 
systems that axe initially prepared in eigenstates of the unperturbed Hamiltonian, by 
assuming that the latter has a continuous spectrum. 

The exponential decay of the survival probability, familiar, for example, in radioac- 
tive decay phenomena, is representative of a purely probabilistic character of the system 
under consideration and is naturally expected to lead to a master equation. This behav- 
ior, however, can be found only at intermediate times, for deviations from it exist both 
at short and long times and can have significant eonsequenees. 

After a short introduction to the long history of the research on the temporal behav- 
ior of such quantum mechanical systems, the short-time behavior and its controversial 
consequences when it is combined with von Neumann's projection postulate in quantum 
measurement theory are critically overvicwcd from a dynamical point of view. We also 
discuss the so-called quantum Zeno effect from this standpoint. 

The behavior of the survival amplitude is then scrutinized by investigating the ana- 
lytic properties of its Fourier and Laplace transforms. The analytic property that there 
is no singularity except a branch cut running along the real energy axis in the first Rie- 
mannian sheet is an important reflection of the time-reversal invariance of the dynamics 
governing the whole process. It is shown that the exponential behavior is due to the 
presence of a simple pole in the second Riemannian sheet, while the contribution of the 
branch point yields a power behavior for the amplitude. The exponential decay form is 
cancelled at short times and dominated at very long times by the branch-point contri- 
butions, which give a Gaussian behavior for the former and a power behavior for the 
latter. 

In order to realize the exponential law in quantum theory, it is essential to take into 
account a certain kind of macroscopic nature of the total system, since the exponential 

behavior is regarded as a manifestation of a complete loss of coherence of the quantum 
subsystem under consideration. In this respect, a few attempts at extracting the expo- 
nential decay form on the basis of quantum theory, aiming at the master equation, are 
briefly reviewed, including van Hove's pioneering work and his well-known "A^T" limit. 

In the attempt to further clarify the mechanism of the appearance of a purely prob- 
abilistic behavior without resort to any approximation, a solvable dynamical model is 
presented and extensively studied. The model describes an ultrarelativistic particle inter- 
acting with A'' two-level systems (called "spins") and is shown to exhibit an exponential 
behavior at all times in the weak-coupling, macroscopic limit. Furthermore, it is shown 
that the model can even reproduce the short-time Gaussian behavior followed by the 
exponential law when an appropriate initial state is chosen. The analysis is exact and no 
approximation is involved. An interpretation for the change of the temporal behavior in 
quantum systems is drawn from the results obtained. Some implications for the quantum 
measurement problem are also discussed, in particular in connection with dissipation. 

1. Introduction and summary 

We know that the radioactive decay is well subject to the exponential law, which 

can be easily understood on the basis of a purely probabilistic argument. We are 
also familiar with this kind of exponential behavior in a dissipative classical system, 
like for example a one-particle motion characterized by friction in an atmosphere 
having a huge number of degrees of freedom. However, an unstable nucleus in 
a radioactive material is usually in an excited quantum mechanical state, whose 
behavior is governed not by classical dynamics but by quantum mechanics. In this 
paper, we review and examine what kind of temporal behavior should take place, 
and up to what extent we can observe the exponential decay in an unstable quantum 
mechanical system. 
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Let US first consider a classical system composed of many particles, in which we 
are interested in a single- mode behavior (for example, a one-particle motion). On a 
fine time scale, whose unit is much smaller than the collision time, we should observe 
a purely dynamical motion governed by the many-body Newton equations. On the 
other hand, if we observe the phenomenon on a macroscopic time scale, whose 
unit is much longer than the typical collision time, we can usually decompose the 
whole motion into two parts, the first being a dissipative systematic motion, and 
the second a fluctuating one, driven by a random force with a vanishing local time 
average. Here we are considering the local time average, defined by an averaging 
procedure over a local time region which is larger than the microscopic relaxation 
time and smaller than the macroscopic time constant. 

In this view, we are dealing with the original many-body problem as a phe- 
nomcnological stochastic process of a single-mode variable (for example, a one- 
particle motion), which is described by the so-called Langevin equation. The sys- 
tematic part, i.e. the local time average of the whole motion, is usually subject 
to the exponential law, which is characterized by a dissipation constant (like the 
friction constant in a one-particle motion). We can derive this description from the 
many-body Newton equations with the help of a coarse-graining procedure. When 
the number of dynamical degrees of freedom of the atmosphere is infinitely large, 
we can obtain the Gaussian distribution of the random force through the central 
limit theorem, and then we can replace, by means of a sort of ergodic theorem^ the 
local time average with a statistical ensemble average for the thermal equilibrium 
with a certain temperature T . 

On the other hand, when the microscopic relaxation time is negligibly small, we 
can assume that the random force is white and then the stochastic process becomes 
Markoffian. We usually call this kind of process a Gauss-Markoffian one. In this 
case, we can easily derive the Fokker-Planck equation, which justifies the approach 
to thermal equilibrium. Usually, the Langevin equation is to be formulated for 
Gauss-Markoffian processes. In this case, the ensemble average of every physical 
quantity can be written as a sum of products of two-time correlation functions. 
Note that we can also ascribe the process driven by a colored noise to a Gauss- 
Markoffian process, if the colored noise can be derived from a white noise through 
a simple mathematical transformation. 

The situation with quantum systems is somewhat delicate for many reasons. In 
order to understand this situation more clearly, let us start our theoretical consid- 
erations from a general discussion on the temporal behavior of quantum-mechanical 
systems. First imagine, for example, that we have a particle system coupled with a 
field in free space. The case of a many-body system with a finite (but very large) 
number of degrees of freedom will be discussed later. 

Many years ago, Weisskopf and WignerS first gave a simple formulation of the 
quantum theory of unstable systems, leading to the famous exponential law. Breit 
and Wigneia introduced an expression for the quantum-mechanical wave functions 
or S'-matrix elements: Such an expression, which is analytic with respect to the 
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energy variable conjugated to time, is equivalent to the exponential law. Even 
though these formulas can work well for phenomenological purposes, the theoretical 
procedure they gave was unfortunately not rigorous. Another important analysis 
yielding the exponential law was given by Gamow in his seminal paper on the 
quantum tunnelling problem. H 

From the point of view of the relationship betweejn the decay phenomena and the 
energy-time uncertainty relation, Fock and Krylovu claimed that the exponential 
decay could not be theoretically accepted. Later, KhalfinEl confirmed this arguments 
on the basis of a mathematical theorem given by Paley and Wiener 

In 1953 HellundS strongly claimed that the long-time behavior is not subject 
to the exponential law, biit rather to a power law. By improving this conclusion, 
Namiki and Mugibayashilj showed that the quantum decay should be subject to 
the following three-step behavior: the Gaussian law at short times, the exponential 
law at intermediate times, and finally the power law at longer tirnes. Later Araki, 
Munakata, Kawaguchi and Goto,§ by making use of the Lee modeljl3 remarked that 
the exponential decay stems from simple poles located on the second Riemannian 
sheet of the analytic expression of the relevant S'-matrix element. A similar but more 
thorough argument was also given by Schwinger.U The quanturn_decay problem was 
also considered and reviewed by Fonda Ghi rardi and RiminiJ13 and was recently 
discussed by Cho, Kasari and Yamaguchill3llj by nmking use of a new solvable model 
of particle decay. Chiu, Sudarshan and Bhamathlla and Horwitzll3 presented other 
solvable models by extending the Lee model. One of the points of great interest in 
Refs. ^ and |l^ was the field-theoretical renormalization procedure. Throughout the 
present paper, however, we shall not enter into the renormalization problem. 

As for a general overview of the analytical expressions, we can refer to early 
work on the dampi ng theory initiated by Heitlerilll See, for example, the paper by 
Arnous and Zienau.Ej One should remember that the time-reversal invariance of the 
whole system is reflected in the fact that the 5-matrix element has no singularity 
on the first Riemannian sheet, except the branch cut running along the real energy 
axis. 

The above-mentioned temporal behavior of unstable quantum systems is closely 
related to the so-called quantum Zeno paradox or quantum Zeno effect in the quan- 
tum measurement problem. The quantum Zeno paradox (QZP), named after the 
famous Greek philosopher Zeno, states that an unstable quantum system becomes 
stable (i.e. never decays) in the limit of infinitely frequent measurements. Of course, 
in practice, we cannot observe this very limit,'' but can only investigate the quantum 
Zeno effect (QZE), i.e. a milder version of the QZP, stating that the probability of 
finding the initial state is increased by a (finite) number of repetitions of a measure- 
ment. We are now simply considering the notion of quantum measurement in terms 
of the naive notion ofjvave function collapse (WFC), namely a simple projection 



onto the initial state.Ea 

The seminal idea of QZP or QZE was introduced under the assumption that the 

^Evcn in principle, this limit cannot be observed, because of the uncertainty principle. See Ref. |l^. 
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Gaussian short-time behavior can be observed and utiHzed in the relevant mi a.n t um 
decay and only the naive WFC takes place in quantum mcasurementsH3^c3 In 
this context, one might think that the observation of this kind of phenomenon 
is a clear-cut experimental evidence in support of the naive WFC. Usually, the 
experimental observation of the Gaussian short-time behavior in a decay process 
is so difficult to perform that CookE3 proposed to use^tomic transitions of the 
oscillatory type, and inspired an important experimented and then an interesting 
debater3't£l As was discussed in the latter papers, we are led to the conclusion 
that the experimental observation of the QZE does not necessarily support the 
naive WFC. In order to explain the situation in detail, we have to examine one 
of the central questions in the quantum measurement problem: What is the wave- 
function collapse? In this context, the issue of the temporal development of quantum 
systems, by which we have obtained the idea of QZP or QZE, is closely related 
to the quantum measurement problem. Two of the present authors have already 
formulated a r e a son able theory of quantum measurements without resorting to the 
naive WFCH'ES In this paper, however, we shall not enter into this kind of problem. 

Let us now consider and discuss a dynamical system composed of a finite number 
of particles put in a finite box. The energy spectrum of this many-body system is 
discrete, and correspondingly, the elastic collision ^-matrix element has a series of 
simple poles running from the minimum energy (say, Eg) to infinity along the real 
axis on the complex E plane. As the number of degrees of freedom and the box size 
become very (infinitely) large in an appropriate way, the energy spectrum becomes 
asymptotically continuous and the series of simple poles approaches a branch cut 
running from Eg to infinity along the real axis, just like in the field-theoretical 
case mentioned above. This is a result of the asymptotic procedure, supplemented 
by an appropriate coarse-graining procedure, that endowes the off-diagonal matrix 
elements of the interaction Hamiltonian with random phases and then provokes the 
appearance of the so-called diagonal singularity^^ A similar procedure was used in 
the theory of nuclear reactions. O As a consequence, the temporal behavior of the 
system is shown to be classified in three steps, the first being Gaussian, the second 
exponential, the third power-like. Usually, the Gaussian period is very short and 
the power decay is very weak, so that the exponential decay seems to dominate over 
the whole process. If we perform a sort of tim&^scale transformation in the weak- 
coupling case, following van Hove's procedure,E3 only the exponential decay part 
remains in the whole process. Being left with the exponential decay only, we can 
regard the whole process as a probabilistic one. Within the theoretical scheme given 
by van Hove, in fact, we can derive the master equation, which is a manifestation 
of irreversible processes, from the quantum-mechanical Schrodinger equation.t^ 

At the beginning of this introductory section, we gave an outline of the deriva- 
tion of the so-called Langevin equation from the many-body Newton equation in a 
classical system. Even in a quantum system, we can follow a similar kind of proce- 
dure by deriving a quantum Langevin equation from the basic quantum-mechanical 
Schrodinger or Heisenberg equations. Strictly speaking, however, we have to notice 
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that aU observables in quantum mechanics are represented by operators, and their 
temporal evolution should obey an operator-valued Heisenberg equation as a first 
principle. Therefore, we should first naturally derive an operator-valued Langcvin 
equation. When and only when the operator nature can be neglected, we can talk 
about a c-number Langevin equation as an approximation. Even in this case, the 
random force never rigorously satisfies the Gauss-Markoffian property, so that we 
cannot calculate all the expectation values of the physical quantities by means of 
a sum of products of two-point correlation functions. When and only when the 
colored property of the quantum random force, and consequently all the important 
quantum properties, can be neglected, we can use a c-number quantum Langevin 
equation, which can be compared with the classical one. 

If we want to formulate a quantum Langevin equation rigorously, we also have 
to take into account the quantum properties of spin statistics, i.e. the well-known 
distinction between fermions and bosons. In order to do this, we should first for- 
mulate an operator-valued Langevin equation for quantum fields, as was already 
proposed by Mizutani, Muroya and Namiki.L3 More details will be published in a 
forthcoming paper. 

It is true that actual physical systems are so complicated that we cannot easily 
obtain exact solutions of the fundamental equation. In this context, it is meaningful 
to examine the temporal evolution of quantum systems by means of exactly sobzahl£ 
models. In this paper, we shall describe in detail some attempts of this kind.E3~Ea 
This is one of the important purposes of the present paper. 

We organize this paper in the following way. In Section 2, the short-time behav- 
ior of quantum mechanical systems is derived, and its relations with the quantum 
Zeno paradox and quantum Zeno effect are critically discussed. After a short in- 
troduction to the classical exponential law and the quantum mechanical deviations 
from it at short times, the difference between QZE and QZP is explained in Sec. 2.1. 
The seminal formulation of the paradox is reviewed in Sec. 2.2, and an explicit ex- 
ample is considered in Sec. 2.3, where it is also shown that both QZE and QZP 
can be given a purely dynamical explanation. It is then shown in Sec. 2.4 that 
the QZP is only a mathematical idealization and is physically unattainable; how- 
ever, its milder version, i.e. the quantum Zeno effect (QZE), can be observed in 
practice. Section 3 deals with the long-time behavior of quantum systems. There, 
the general expression for the survival probability amplitude is given for quantum 
systems with a huge (ideally infinite) number of degrees of freedom, where the so- 
called random-phase approximation seems most plausible. We see that the Fourier 
or Laplace transforms of the amplitude can clarify its analytical structure, which in 
turn determines its temporal behavior. The following subsections aim at giving a 
deeper insight into the profound connection between the analytic structure and the 
time dependence of the amplitude. We also confirm the general statements on the 
basis of a concrete, but still rather general Hamiltonian. After a naive (textbook) 
derivation of the the short-time Gaussian behavior and the long-time exponential 
law, a counter theorem claiming the unattainability of the exponential form at very 
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long times is introduced in Sec. 3.1, and the analytical structure of the survival 
amplitude is reinvestigated carefully and explicitly in Sec. 3.2. It is shown that 
the exponential behavior is due to a simple pole located on the second Riemannian 
sheet, and turns out to be dominated by a power decay at very long times. These 
exponential and power-like behaviors are a manifestation of the analytic property 
of the amplitude and the role played by the branch cut in the complex energy plane 
is further analyzed in Sec. 3.3 on the basis of a concrete Hamiltonian. Section 4 is 
devoted to an explicit demonstration of the exponential decay form on the basis of 
a solvable dynamical model. The model is solved exactly and a propagator, repre- 
sentative of the survival amplitude in this case, shows the exponential form at all 
times in the weak-coupling macroscopic limit. The temporal behavior is shown to 
be sensitive to the choice of the initial state and the short-time Gassian as well as 
the following exponential behaviors are derived in the case of a wave-packet initial 
state. A possible relation between the notions of dissipation and decoherence is 
mentioned, in relation to the quantum measurement problem. The last section. 
Sec. 5, is devoted to conclusions and additional comments. 

2. Short time behavior and quantum Zeno effect 

We shall start by giving a brief outline of the problem. In classical physics, an 
expression for the decay probability of an unstable system is easily obtained by a 
heuristic approach (see, for example, Ref. |^): One assumes that there is a decay 
probability per unit time T that the system will decay according to a certain specific 
process. Such a probability per unit time is constant, and does not depend, for 
instance, on the total number N of unstable systems, on their past history, or on 
the environment surrounding them. If the number of systems at time t is N(t) (a 
very large number), the number of systems that will decay in the time interval dt 
is 

dN 

dN = -NTdt or = -TN, (1) 

dt 

which yields 

Nit) = 7Voe~", (2) 
where A^o = A^(O)- One defines then the "survival" or "nondecay probability" 

PW^^^e- (3) 

The (positive) quantity T is interpreted as the inverse of the "lifetime" r. Notice 
the short-time expansion 

P{t) = 1-Tt + ---. (4) 

Even though the above derivation can be found in elementary textbooks at un- 
dergraduate level, the assumptions underpinning it are delicate. The reader will 
recognize the essential basic features of a Markoffian process, in which memory 
effects are absent. More to this, one is excluding a priori the possibility that co- 
operative effects take place, making F and P environment-dependent. It is not 
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surprising, then, that the resultant solution (g) exhibits all the basic ingredients of 
a dissipative behavior. 

Let us now turn to a quantum mechanical description of the phenomenon. Let 
|a) be the wave function of a quantum system Q at time t = 0. The evolution of Q is 
governed by the unitary operator U{t) — exjp{~iHt), where H is the Hamiltonian. 
We define the survival or nondecay probability at time t as the square modulus of 
the survival amplitude 

P{t) = \{a\e-^"'\a)\\ (5) 
A naive and elementary expansion at short times yields 

Pit) = l-tma\H^\a)^i(a\H\a)r)+--. 

= l-t^AHf + ---, (6) 

which is quadratic in t and therefore yields a vanishing decay rate for i — > 0. Observe 
that we are implicitly assuming that \a) be normalizable and (Ai?) nonvanishing, 
or, in other words, that \a) is not an eigenstate of H. (Note that one simply gets 
P{t) — 1 if |a) is an eigenstate of H.) We also assumed that all moments of H in the 
state \a) be finite. Under these assumptions, we can easily infer that the survival 
probability at short times is of the Gaussian type: 



P(t)c.exp(-- : T^'^iAHf. (7) 



A 

The temporal behavior of P{t) is closely related to the famous issue of the time- 
energy uncertainty relation, which has been a long-standing argument of discussion. 
In this context, we can refer to old papers, like for instance those by Fock and 
Krylovi and Khalfin.0 The former claimed that the exponential decay was not 
theoretically acceptable, and the latter confirmed this conclusion on the basis of a 
mathematical theorem given by Paley and Wiener .i 

Needless to say, the above result is in manifest contradiction with the exponential 
law (H) that predicts a constant, nonvanishing decay rate F. More to this, there are 
some other important differences between the classical and the quantum mechanical 
results. First of all, observe that the derivation of holds for individual 

systems, while that of (^) makes use of an ensemble made up of a huge number 
N of identically prepared systems. This is because probabilities are ontological 
in quantum mechanics, unlike in classical physics. Second, in the quantum case 
the initial state plays an important role and directly enters in the definition of P, 
raising subtle questions about the problem of the initial state preparation, as is 
widely known: See, for example, Refs. |3^ and ^ By contrast, in the classical case, 
the initial state plays no crucial role, and simply characterizes the constant A^o in 
(||). This is obviously related to the underlying MarkofRan assumptions. Finally, 
notice that the quantum mechanical analysis holds true for any initial state |a), not 
necessarily unstable. 

The requirements that the initial state \a) be normalizable and its energy finite 
are very important: If these two conditions are not met, the behavior of the survival 
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probability can be different from (p) or (M). A typical example is the Breit-Wigner 
spectrum, as was discussed in Refs. p| and |39|. The above two requirements will play 
an important role in the model we shall study in Sec. ^. 

We shall postpone our considerations on the temporal behavior at longer times 
to Sec. ||. Now we analyze one of the most striking and famous consequences of the 
quadratic short-time behavior: The quantum Zeno effect. 

2.1. Quantum Zeno paradox and quantum Zeno effect 

Zeno was born in Elea, a southern Italian town not far from Naples, about 490 
B.C. At those times, that region was part of the so-called Magna Graecia, and was 
culturally Greek. Zeno was a disciple of Parmenides, the most prominent figure 
of the so-called Eleatic school of philosophers. Parmenides firmly believed in a 
unique, entire, total Truth ("being"), and could not accept the idea that this Truth 
could change ("becoming") or be composed of smaller entities ("multiplicity"). His 
dialogue with Socrates on these matters is one of the most famous debates in the 
history of Greek philosophyJlj 

Even though Zeno's contribution to the philosophy of "being" were not as impor- 
tant as Parmenides's, the former was an unequaled orator, and is believed to have 
invented dialectic, that was subsequently so widely used by Socrates. Two argu- 
ments given by Zeno in support of Parmenides's philosophical viewpoints are most 
famous. In the first one, Achilles cannot reach a turtle because when the former 
has arrived at the position previously occupied by the latter, the turtle has already 
moved away from it, and so on ad infinitum. In the second one, that concerns us 
more directly, a sped arrow never reaches its target, because at every instant of 
time, by looking at the arrow, we clearly see that it occupies a definite position in 
space. At every moment the arrow is therefore immobile, and by summing up so 
many "immobilities" it is clearly impossible, according to Zeno, to obtain motion. 
If he lived nowadays, Zeno would probably consider the existence of photographs 
(in which all moving objects are still) as the best proof in support of his ideas. 

In the light of his paradoxical arguments, Zeno can be considered a forerun- 
ner of sophism. We know today that the two above-mentioned paradoxes can be 
solved by infinitesimal calculus. Nevertheless, one should not overlook the fact that 
Zeno basically aimed at giving very provocative arguments against the concept of 
"becoming" , in order to ridicule the critics who tried to deride the philosophy of 
"being" . 

It is somewhat astonishing that a conclusion similar to Zeno's holds in quantum 
theory: It is indeed possible to exploit the quantum mechanical vanishing decay 
rate at short times (^) in a very interesting way, by slowing down (and eventually 
halting) the deca^i^ process. The essential features of this effect were alreadi!_known 
to von Neumanrll^l and were investigated by several authors in the past.E3~E3 Misra 
and Sudarshan first named this phenomenon after the Greek philosopher. 

We shall first give an elementary derivation of this effect, and present a more 
rigorous result in the following subsection. Let us start from Eq. (^) or (0), that we 
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rewrite as 

P(i) = l-4 + ---=^e-*V-S (8) 

where tq is the characteristic time of the Gaussian evolution. Suppose we perform 
N measurements at equal time intervals t, in order to ascertain whether the system 
is still in its initial state. After each measurement, the system is "projected" onto 
the quantum mechanical state representing the result of the measurement, and the 
evolution starts anew. The total duration of the experiment is T = Nt. The 
probability of observing the initial state at time T, after having performed the N 
above-mentioned measurements, reads 

Notice that both T and N are finite, in the above. This is the quantum Zeno effect: 
Repeated observations "slow down" the evolution and increase the probability that 
the system is still in the initial state at time T. In the limit of continuous observation 

{N oo) one obtains the quantum Zeno paradox 

pW(T)^h--^(|) j "^""l- (10) 

Infinitely frequent observations halt the evolution, and completely "freeze" the ini- 
tial state of the quantum system. Zeno seems to be right, after all: The quantum 
"arrow" , although sped under the action of the Hamiltonian, does not move, if it is 
continuously observed. 

What is the cause of this peculiar conclusion? Mathematically, the above result 
is ascribable to the general property 

pi^)iT)^(i-o^j,^y ''-^i. (11) 

However, the above limit is unphysical, for several reasons. Indeed, notice that there 
is a profound difference between the A^-finite and the A^-infinite cases: In order to 
perform an experiment with N finite one must only overcome practical problems, 
from the physical point of view. (Of course, this can be a very difficult task, in 
practice.) On the other hand, the N oo case is physically unattainable^ and is 
rather to be regarded as a mathematical limit (although a very interesting one). 
In this sense, we shall say that the quantum Zeno effect, with A^ finite, becomes a 
c^uantum Zeno paradox when A^ — s- cx).E3 

It must also be emphasized that in the above analysis the observations (measure- 
ments) are instantaneous, and this is why it is possible to consider the N oo limit: 
No time is spent in measurements.EEl This is in line with the Copenhagen school of 
thought, and can be regarded as a common and general feature of von Neumann- 
like descriptions of a measurement process. Even though such a picture is often 
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accepted among physicists, it is very misleading, in our opinion, and hides some 
very important aspects of the problem: Indeed a measm'emcnt process, if analyzed 
(as it should) as a concrete physical process, takes place during a very long time on 
a microscopic scale, although we can regard it as if it happened instantaneously on 
a macroscopic scale. 

The physical unrealizability of the above-mentioned limit will be thoroughly 
discussed in Sec. where it will be shown that, as a consequence of the uncertainty 
relations, the N oo limit turns out to be impossible, even in principled 

2.2. Misra and Sudarshan's formulation of the quantum, Zeno paradox 

Let us now look at the seminal formulation of the quantum Zeno paradox as given 
by Misra and Sudarshan.@ The derivation of the paradox is entirely based on von 
Neumann's projection rule and therefore hinges upon the concept of quantum mea- 
surement process. The reader should notice the fundamental role played by the 
projection operator O in the following, and should remember that its effect on the 
state of the system is supposed to take place instantaneously. We shall argue, in 
the next subsection, that the quantum Zeno effect and paradox can be derived on 
a purely dynamical basis, without making use of von Neumann's projections .E3 

Q is our unstable quantum system, whose states belong to the Hilbert space 
Ti. and whose evolution is described by the unitary operator U{t) — exp{~iHt), 
where H is a semi-bounded Hamiltonian. The initial density matrix of system Q 
is assumed to be an undecayed state po, and let O be the projection operator over 
the subspace of the undecayed states. By definition, 

Po^OpoO, Tr[poO] = l. (12) 

Assume that we perform a measurement at time t, denoted by the projection oper- 
ator O, in order to check whether Q is still undecayed. Accordingly, the state of Q 
changes into 

Po^pit)^OU{t)poUHt)0, (13) 
so that the probability of finding the system undecayed is given by 

p{t)=TY[U{t)poU^{t)0]. (14) 

The process (|l^) will be referred to as "naive wave function collapse" . As already 
stressed, the observations (measurements) schematized via the operator O take 
place instantaneously. 

The formulation of the QZP proceeds as follows: We prepare Q in the ini- 
tial state Po at time (this is formally accomplished by performing an initial, 
"preparatory" measurement of O) and perform a series of observations at times 
T/N, 2T/N, ...,{N- 1)T/N, T. The state p'-^^T) of Q after the preparation and 
the above-mentioned N measurements reads 



pW(r) = lV(r)pol4(r), VN{T)=[OU{T/N)Of (15) 
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and the probability of finding the system undecayed is given by 

pW(0=Tr[y^(T)po^^(T)] . (f6) 

So far, the derivation is straightforward^ and does not involve, in our opinion, any 
new concept. By contrast, much more delicate is the set of assumptions leading 
to the "paradox" : One considers the so-called limit of continuous observation, in 
which iV — !■ oo, and defines 

V(T) EE hm V^iT), (17) 

provided the above limit exists in the strong sense. The final state and the proba- 
bility of observing the undecayed system read then 

p{T) = V(T)poVt(r), (18) 
V{T) EE lim P(^)(r) ==TrfV(T)/9oV^(r)l . (19) 

Moreover one assumes, on physical grounds, the strong continuity of V{t) in i = 0, 

lim V(t) = O. (20) 

Misra and Sudarshan then proved that under general conditions'^ the operators 
V(T) (exist for all real T and) form a semigroup labeled by the time parameter T. 
Moreover, Vl'(r) = V(-r), so that Vl'(r)V(T) ^ O. This implies, by virtue of Eq. 
(H), that 

V{T) ^ lim P(^)(r) = Tr [poO] = 1. (21) 

If the particle is continuously observed (to check whether it decays or not), it is 
"frozen" in its initial state and will never be found to decay! This is the essence of 
the "quantum Zeno paradox" . 

Two important problems are still open, at this point. First: Is the present 
technique, that strongly hinges upon the action of projection operators d la von 
Neumann, really necessary in order to derive the QZP or the QZE? Second: Is the 
N ^ oo limit physically sensible? 

We shall tackle the first question in Sec. where a purely dynamical derivation 
of the QZP and QZE is presented, that makes use of unitary evolutions and does 
not involve any projection (provided a final measurement is carried out). We shall 
come back to the second question in Sec. where the N ^ oo limit will be shown 
to be in contradiction with the Heisenberg uncertainty principle. 



Except, of course, that one is led to wonder about the physical meaning of the operator O, 
and about the reason why a "projection", unlike all other physical phenomena, should occur 
instantaneously. These questions are left to the reader. The von Neumann formalism was critically 
analyzed in Refs. Bll and B2l 

'^Somc technical requirements are needed for the proof. 
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2.3. Dynamical quantum Zeno effect: A simple model 

Let us show that it is possible to otain both the QZE or jl^ ) and the QZP ( ]To| ) 
or (|2^) by making use of a purely dynamical process. The general case is treated 
in Ref. |3^. Here, we shall not give a general proof of the above statement, because 
the best proof is by inspection of an explicit example: Consider the experimental 
setup sketched in Figure 1. An incident neutron, travelling along the y-direction. 



Fig. 1. (a) "Free" evolution of the neutron spin under the action of a magnetic field. An emitter E 
sends a spin-up neutron through several regions where a magnetic field B is present. The detector 
Do detects a spin-down neutron: No Zeno effect occurs, (b) The neutron spin is "monitored" at 
every step, by selecting and detecting the spin-down component. Dq detects a spin-up neutron: 
The Zeno effect takes place. 



interacts with several identical regions in which there is a static magnetic field B, 
oriented along the x-direction. We neglect unnecessary complications and describe 
the interaction by the Hamiltonian H = fiBai (fi being the (modulus of the) neutron 
magnetic moment, and cti the first Pauli matrix.) Let the initial neutron state be 
Po = PtT = I t)(T 1 1 where | t) and | J,) are the spin states of the neutron along 
the z axis, which can be identified with the undecayed and decayed states of the 
previous subsections, respectively. Assume that there are N regions in which B is 
present and that the interaction between the neutron and the magnetic fields has 
a total duration T = Nt [t = £/v, where i is the length of each region where B is 
present and v the neutron speed). It is then straightforwardEj to show that if T is 
chosen so as to satisfy the "matching" condition T = (2m + 1)tt /uj, where m is an 
integer and uo = 2fj,B, the final state and the probability that the neutron spin is 
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found down at time T read respectively 

PiT) = Pu, (22) 
Px(T) = 1. (23) 

The experimental setup in Figure 1(a) is such that if the system is initially prepared 
in the up state, it will evolve to the down state after time T (7r-pulse). 

The situation outlined in Figure 1(b) is very different. The experiment has been 
modified by inserting at every step a device able to select and detect the down 
component of the neutron spin. This is accomplished by a magnetic mirror M 
and a detector D. The magnetic mirror yields a spectral decompositio niES by 
splitting a neutron wave with indefinite spin (a superposed state of up and down 
spins) into two branch waves (each of which is in a definite spin state along the z 
axis) and then forwarding the down component to a detector. The action of the 
magnetic mirror can be compared to the inhomogeneous magnetic field in a typical 
Stern-Gerlach experiment. It is very important, in connection with the QZE, to 
bear in mind that the magnetic mirror does not destroy the coherence between the 
two branch waves: Indeed, the two branch waves corresponding to different spin 
states can be split coherently and brought back to interfere.c^ The global action of 
M and D can be formally represented by the operator O = p^p If the initial Q 
state and the "matching" condition for T = Nt are the same as before, the density 
matrix and the probability that the neutron spin is up at time T — (2m + 1)tt/uj 
read respectively 

pW(T) = VMiT)poVliT)^^cos'^j p^^^(cos2^) (24) 

This discloses the occurrence of a QZE: Indeed, p|^^(r) > P^'^^^\t) for > 2, 
so that the evolution is "slowed down" as N increases. In the limit of infinitely 
many observations 

^W(r) piT)^Pn, (26) 

pj^)(r) 7't(T) = 1. (27) 

Frequent observations "freeze" the neutron spin in its initial state, by inhibiting 
{N > 2) and eventually hindering (TV oo) transitions to other states. Compare 
Eqs. ( ^6|) and ( ^ ) with ( p2|) and (|2^): The situation is completely reversed. 

It must be observed, however, that it is possible to obtain the same result with- 
out making use of projection operators, by simply performing a different analysis 
involving only unitary processes. Observe first that the preceding analysis involves 
only the Q states. If the state of the total (neutron + detectors) system is duely 
taken into account, the final total density matrix, in the channel representation. 
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reads 



/ 



0\ 



V 







i,j = 0,l,...,iV, (28) 



where c — cos(7r/2iV) and s = sin(7r/2A^). This corresponds to the case of frequent 
observations, in which the neutron route was observed at every step. The i = j = Q 
component corresponds to detection by Z3oj while the i = j ^ n (n — 1, . . . , N) 
component corresponds to detection in channel N — n + 1. Observe that in the 
above expression the total density matrix has no off-diagonal components as a 
consequence of the "wave function collapse" by measurement. 

Remove now Di, • • • , Dat in Figure 1(b): In other words, no observation of the 
neutron route is carried out (except the final one performed by Dq). A straightfor- 
ward calculation yields the following final density matrix for the Q system 



V 



r,2N 



.2N~1 ■J2N~2 



-ISC 



„2N-2 



ISC 

s^c^'^-^ s^c^'^-^ 



ISC 



N 



hj = 0, 1, 



J 



The two expressions above clearly show that we have the same probability p|^'' = 
[cos2(7r/2JV)]''^ of detecting a spin-up neutron at Dq irrespectively of the presence 
of detectors Di, . . . , Dn in Figure 1(b). It appears therefore that no projection 
rule is necessary in this context. The quantum Zeno effect can be given a purely 
dynamical explanation. 

If iV — > oo, both density matrices tend to the limiting expression 



700 



/ 1 





V ; ; : 



i,j = 0,1,. 



(30) 



This is the quantum Zeno paradox. It can be obtained both by making use of 
projection operators a la von Neumann and by means of a purely dynamical pro- 
cess. The above conclusions can be generalized to an arbitrary quantum system 
undergoing a Zeno- type dynamics, and it can even be shown that it is possible to 
"mimic" the instantaneous action of a projection operator by making use of the 
impulse approximation in quantum mechanicsEj 

It must also be stressed that an idea analogous to the one described in this 
subsection was outlined by Peres0 who made use of photons, rather than neutrons. 
His proposal inspired an interesting experiment O 
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2.4. The N ^ oo limit and its physical unrealizability 

Let us now discuss the meaning of the iV — s- oo hmit and show that it is unphysical0 
We start by observing that the condition uiT = LoNt = (2TO+l)7r, which is to be met 
at every step in Figure 1(b), imphes (by setting m = 1 for simphcity and without 
loss of generahty) 

where all quantities were defined before Eq. (^). Obviously, as N increases in the 
above equation, the practical realization of the experiment becomes increasingly 
difhcult. But close scrutiny of Eq. ( ^ ) shows that p'^^\t) cannot tend to 1, even 
in principle, in the N oo limit, because of the uncertainty relations. Indeed, let 
(j) be the argument of the cosine in Eq. (p5|): 



(32) 



Mathematically, the above quantity is of order N^^. On the other hand, from a 
physical point of view, it is impossible to avoid uncertainties in the neutron position 
Ax and speed Av. As a consequence, (f) cannot vanish, because it is lower bounded 
as follows 

hvQ hvo 2MvqAv^ 
where M is the neutron mass and we assumed that the size i of the interaction region 
(where the neutron spin undergoes a rotation under the action of the magnetic field) 
is large r th an the longitudinal spread Ax of the neutron wave packet. An accurate 
analysialj shows that the same bound holds in the opposite situation (£ < Ax) as 
well. 

By defining the magnetic energy gap AE„i = 2^B and the kinetic energy spread 
of the neutron beam AEk ~ A{Mv'^ /2)\y=yg, the above inequality reads 



1 AE„ 
> -- 



4 AEk 



(34) 



It is now straightforward to obtain an expression for the value of the probability 
that a spin-up neutron is observed at Dq when N is large: 



P{^)(T) 



(cos 00 ) 



2JV 



2N 



1 



1 f AEr, 



32 V AEk 



2N 



(35) 



Notice that not only the above quantity does not tend to 1 , but it vanishes in the 
N oo limit. In other words, in the experiment outlined in Figure 1(b), no spin-up 
neutron would be observed at Dq in the N oo limit! 

What is it reasonable to expect, in practice? Even though the analysis of the pre- 
vious subsections does not take into account the limits imposed by the uncertainty 
principle, it must be considered that, in practice, N cannot be made arbitrarily 
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large. In order to evaluate how big N can be in order that the QZE be observable 
in the above experiment, set P^^\t) ^ 1/2. We get 

641n2 , , , 

where we assumed reasonable values for the energies of a thermal neutron. In 
conclusion, N turns out to be large enough in order that the QZE be experimentally 
observable, at least up to a certain approximation. 

Criticisms against the physical meaning and realizability of the N oo limit 
were put forward some years ago by Ghirardi et ahEl Although different from ours, 
these criticisms were based on the time-energy uncertainty relations. Our argument, 
outlined for neutrons, holds true in general, and even in the recent experiment 
performed with photons in Innsbruck.Ea An exhaustive analysis is given in Ref. |l^. 

There are in fact other strong arguments against the iV — > oo limit, from a 
physical point of view. For instance, the above calculation considers only the time 
spent by the neutron in the magnetic field B. In practice, however, one cannot 
neglect the time elapsed during the interaction between the neutron and a magnetic 
mirror M, which is of the order of 10^^ — 10^^ s. 

The above discussion should have made it clear that although the short-time 
behavior is essentially governed by an evolution of the Gaussian type, the limit of 
continuous observation {N oo) should never be lightheartedly taken for granted. 
It is safer and much more reasonable to start one's considerations from formulas 
like (^), where both T and N are finite. 

On the other hand, we shall soon see that there are physical situations in which 
the temporal behavior of a quantum system can be considered approximately expo- 
nential (the mathematical and physical conditions of validity of such an assumption 
will be thoroughly discussed in the following section). In such a case, the survival 
probability at time T after TV measurements would read 

where the subscript e stands for "exponential" , and te is the characteristic time of 
the exponential law. Notice that, unlike in (|^), in this case the final result does not 
strongly depend on N. 

If there exist finite values of T and N such that 

T T"^ T Tr 

TE ~ T^N Ntg~ te' ^ ' 

then the exponential and the Gaussian regions are comparable with each other. 
Which temporal behavior is actually to be observed depends therefore on the rel- 
ative magnitude of the parameters tq and te- For this reason we may not com- 
pletely exclude the possibility of observing a sort of quantum Zeno effect within the 
framework of the exponential decay. An analogous point was raised by Schulman, 
Ranfagni and Mugnaijf] within the context of the WKB approximation. 
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It should be emphasized that it is very difficuh to obtain general estimates of the 
two characteristic times mentioned above. For example, Khalfin, as well as other 
authors, even considered and critically discussed the possibility that the proton 
decay has never been observed because its Gaussian characteristic time t^°^ might 
be longer than the lifetime of the Universe.C3 

3. Quantum behavior at longer times 

In the first part of the preceding section, we have seen that the temporal behavior 
of the survival probability at short times is well described by a Gaussian form. 
Our main task in this section is to investigate the temporal behavior of quantum 
systems at longer times and to understand the status of the familiar exponential 
law in quantum theory. 

Since an articulate and complete analysis of the temporal behavior of quantum 
mechanical systems, especially at very long times, is a rather complex issue, we 
shall organize our discussion in the following way. First, in Sec. |.|l|, we consider a 
particular interaction Hamiltonian and discuss the temporal behavior both at short 
and long times. Then, in Sec. |^.^, we reconsider our conclusions in the most general 
case and discuss which problems arise in the very long time region as a consequence 
of the analytic properties of the survival amplitude. Finally, for the sake of clarity, 
in Sec. we return to the particular case and thoroughly discuss the very long 
time behavior. 

Suppose that the total Hamiltonian is decomposed into two parts, H = Hq + H', 
where Hq and H' arc the unperturbed and the interaction Hamiltonians, respec- 
tively. Let \n) be an eigenstate of the former belonging to the eigenvalue En 

Ho\n)^En\n), l = (1) 

n 

Even though we use the discrete notation for notational simplicity, it is to be un- 
derstood that the possibility of continuous spectrum is included as well. At t = 0, 
the interaction is turned on and the system starts to evolve from an initial state, 
say \a), into the other states. The initial state is not an eigenstate of H, and is 
prepared as an appropriate wave packet, that can be represented by a superposition 
of |n)'s. For the sake of simplicity, we shall denote the initial state by an eigenstate 
of Hq, say \a). The interaction Hamiltonian H' is chosen so that the condition 

{n\H'\n) = (2) 

is satisfied Vrt. This ensures that the Lippmann-Schwinger equation 

l^a) = \a) + ^ H'\4>a) (3) 
£ja — 

can be solved. Notice that the condition corresponds to mass rcnormalization 
in field theories. As already mentioned, we shall start by discussing the temporal 
behavior of the survival probability amplitude in a simple case. 
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3.1. Derivation of temporal behaviors: Standard textbook method 

Here we shall choose a particular interaction Hamiltonian H' subject to the condi- 
tion that it has nonvanishing matrix elements only between the initial state \a) and 
the other states \n) [n ^ a): 

{n\H'\n')^Q, {a\H'\n)^Q, {n,n' ^ a). (4) 

We shall see later that the first relation can be considered as a stronger version of the 
so-called random phase approximation (see (^8|) below). Observe that these condi- 
tions are peculiar to some solvable models. Historically, the Lee model was based on 
a similar idea.EJ Nevertheless, the above can still be considered as a rather general 
form of the interaction Hamiltonian, as can be seen in some standard textbooko 
We arc interested in the temporal behavior of the survival probability amplitude 
{a\Ui{t)\a) . The evolution operator in the interaction picture Ui{t) satisfies the 
Schrodinger equation 

i^Uiit) - H',{t)Ui{t), H'jit) - e'^°*i7'e^'^"*, (5) 
at 

where H'j[t) is the interaction Hamiltonian in the interaction picture. The solution 
Ui{t) can be obtained iteratively from the recursion relation 

Ui{t) = l~i I dt'Hi{t')Ui{t') 
Jq 

= 1 - z / dt'H'jit') + i-tf f dti f dt2H'j{ti)H'j{t2)Ui{t2), (6) 

Jo JQ Jq 

under the initial condition ?7/(0) = 1. Notice that the conditions (||) and (^) for 
the interaction Hamiltonian imply 

{a\H'j{tM{t2)Ui{t2)\a) = {a\H'i{ti)H'j{t2)\a){a\Ui{t2)\a). (7) 

This relation enables us to obtain a closed equation for the diagonal matrix element 
of the evolution operator {a\Ui{t)\a) (survival amplitude) 

{a\Ui{t)\a) = 1- f dt^ r di2e'^"^*^-*^V(ii - t2){a\Ui{t2)\a), (8) 
Jq Jq 

or an integro-differential equation 

j^{a\Ui{t)\a)=- j\t^e'''-^'-'^^f{t-t^){a\Ui{t^)\a) (9) 

with the initial condition {a\Ui{Q)\a) ~ 1. Here the function f{t) is defined as 

f{t) = {a\H'e-^"°'H'\a) . (10) 

Note that (^ is a series of repeated convolutions of e*^°^*i~*=' V(^i — ^2), and all 
other higher-order terms disappear, due to the particular conditions (^. From this 
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point of view, it is convenient to reduce or to an algebraic equation by means 
of a Laplace transform, which is quite appropriate in order to incorporate the initial 
condition. The solution can now be written as an inverse Laplace transform 

1 pico+e s 

{a\Uiit)\a) = — / ^^ds, (11) 
where the function g is given by 

/"OO 

g(s,t) = s + t e-'"e*^""/(w)du. (12) 



(The second term in g(s, t) will turn out to correspond to the self-energy part l^aiE) 
in the general case. See the next subsection.) 

Having obtained the solution in a closed form, we can discuss its temporal 
behavior and extract its explicit time dependence both at short and long times. For 
small t, the function g is expanded around t = 



g(s,t) = s + t^ e-'''e'^'^'''f{tu)du 
Jo 



/(0),2 



s 

which determines the behavior at small t 



du 

i\ (13) 



f(0) 9 /(0),2 ^ 

~ 1 - - e — ^* . (14) 

Notice that /(O) = {a\H''^\a) is positive definite. The system exhibits the Gaussian 
behavior for small t, in accordance with the general argument stated in Sec. |^. 
[Specialize Eq. (^) in Sec. ^ to the case in which \a) is an eigenstate of Hq, under 
the condition (0).] 

On the other hand, for large t, we can see that the amplitude {a\Ui{t)\a) behaves 
quite differently. In this case, the function g is expanded as a power series of 1/t 

g{s,t) = s + t I e-'"e'^""/(M)du 

(l-^u+---) e''^^''f{u)du 
sll- I ue'^''''f{u)du)+t[ e'^'^''f{u)du. (15) 



Therefore we obtain the following exponential behavior at longer times 

{a\Uj{t)\a) ^ e-^\ (16) 
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as long as the exponent A, given by 

exists. The derivation explained here is similar to that given by Messiah,ii even 
though we endeavoured to bridge the gap with the general theory to be sketched 
in the following subsection. Notice that the above form (|l6|) does not necessarily 
mean that the survival probability decays exponentially. In order to obtain such a 
behavior we have to show that the real part of the exponent A is positive. 

The positivity of the real part of the exponent A cannot always be shown, as can 
be seen in some simple solvable models, where only the oscillating behavior appears 
at all times. Consider, for example, the Hamiltonian H = Hq + H' = fiBo(Ts + fiBai, 
where Hq — jiB^aj, is regarded as the free part and H' = fiBai as the interaction 
(a similar example was considered in the context of the quantum Zeno effect, in 
the paper by Inagaki et al.cj). If we start from the initial state \a) = | t); the 
above analysis yields (see ( p^ for the definition of /) f{t) — (/ii3)^e*'^^''* and 
a purely imaginary A = i^BoB^ /{2Bq + B"^). As was to be expected, Eq. ( |l6| ) 
yields an oscillatory behavior at all times and no exponential law can be seen. 
This is ascribable to the finiteness of the system considered: We cannot expect any 
exponential behavior if the number of degrees of freedom of the system is finite. As 



a matter of fact, the numerator appearing in A in (17) is purely imaginary unless 



the initial state \a) is continuosly degenerated w.r.t. energy. We have 

/ e'^-^'fiu^du = lim (a\H' H'\a) 

Jo ^-+0 Ea-Ho+ie 

= {a\H'P—^—H'\a) + ix{a\H'5{Ea - Ho)H'\a) 

■ \^ \{n\H'\a)\^ , ^ , 2 /TON 

= * F _ F +^ \{n\H\a)\. (18) 

This expression implies that if \a) is not degenerated the last term gives no con- 
tribution, which makes the above quantity and therefore the exponent A purely 
imaginary. The degeneracy of the initial state is necessary for the survival probabil- 
ity to decay exponentially. Remember that we are mainly concerned with a particle 
system coupled with fields in free space or with an interacting system made up of 
a particle and a many-body system with a huge number of degrees of freedom, put 
in a finite but very large box, as was mentioned in the introductory section. For 
such systems, we naturally expect the continuous energy spectrum to be always 
degenerated. 

It is worth_noting that the above expression is a realization of the famous Fermi 
Golden Rule,Ey according to which the decay rate is given by 
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by first-order perturbation theory. We understand that twice the real part of A, 
standing for our decay rate, coincides with the above T in the weak-interaction limit 
(so that the denominator in (^) is approximately equal to one). Needless to say, 
this formula can be improved by including higher-order perturbation terms in a 
well-known way. 

Although the condusion we have obtained seems reasonable and its derivation 
sound and standardjlj it is known that at very long times quantum systems never 
show the exponential behavior, a fact in contradiction with the above result as well 
as with our naive expectation. In fact, the unattainability of the exponential law at 
very long times is a mathematically unavoidable consequence. In order to see this, 
recall that the evolution operator in the interaction picture Ui{t) is expressed as 

Ui{t)^e'"°'e"'"\ (20) 

in terms of the total Hamiltonian H . By introducing a complete orthonormal set 
of eigenstates of the latter, we easily arrive at the following form for the survival 
probability amplitude 

{a\Ui{t)\a) = e^^"* j uja{E)e~''''dE, (21) 
where LUa{E) is the energy density of the initial state 

cu^iE)^ J2 \{n\a)\'- (22) 

Here the summation is taken over all the quantum numbers, except energy, that 
are necessary for the specification of a complete orthonormal set. We see that 
the survival probability amplitude is the Fourier transform of the energy density 
uja{E), a fact first pointed out by Fock and Krylov.Q Now let us assume, on physical 
grounds, that the spectrum of the total Hamiltonian is bounded from below so that 
the vacuum state is stable: There exists a, certain finite energy Eg below which the 
function uJa{E) vanishes. Khalfin showedQ that this condition on cua requires that 
its Fourier transform (the survival probability amplitude) satisfies the inequality 

|ln|(a|[/Kt)|a)|| ,, ^ 

as a consequence of the fundamental Paley- Wiener Theorem.! The inequality (^ ) 
implies that the survival probability does not decay exponentially: The decay pro- 
cess proceeds more slowly than exponentially at large times. Notice that the only 
assumption made in the above argument is the existence of a finite Eg: Its very 
value is irrelevant. Therefore the conclusion is quite general.'^ 

At this point, we may feel somewhat puzzled with the above results, especially 
with that at long times, because we are familiar with the exponential decay law in 

'^A full account of this theorem, related topics and further references are given in Ref. 
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classical theory and naturally expect the same form to be valid also in quantum 
theory. Indeed, the exponential law for decay processes has been well confirmed 
experimentally: No deviation from it at long times has ever been observed. A list 
of experimental results, examining the exponential law, can be found in Ref. 
Then how can we reconcile quantum theory with the exponential law? 

The above-mentioned theorem ( |2^ ) requires that the amplitude decays more 
slowly than exponentially in the large-time limit t oo. Therefore we may expect 
it to decay exponentially at intermediate times, even if not at very long times. At 
any rate, something must have been overlooked in our derivation of (p^. In the 
next subsection, we shall reanalyze the temporal behavior in full generality and 
shall focus our attention on the analytic properties of the amplitude: Indeed, in 
the above analysis, we have assumed that the order of the integrations over s and 
u and the i — > and t ^ oo limits can be interchanged, which is licit only when 
the integrals converge uniformly and the initial state is a well-behaved wave packet. 
This is by no means trivial, and is where subtleties might come into play. 

3.2. General arguments 

Let us analyze in full generality the temporal behavior of the survival probability 
amplitude {a\Ui(t)\a) . The evolution operator in the interaction picture Ui{t) is 
related to its counterpart U{t) in the Heisenberg picture and to the S'-matrix in the 
usual way: 

U{t) = e-'"* = e-'"°'Ui{t) e-'"°'S , (24) 

where H and _ffo stand for the total and the free Hamiltonians, respectively. It 
is important to notice that the unitary operator U{t) is expressed as a Fourier 
transform 

U{t) I e-^^'G{E)dE; G{E) = — L- , (25) 

and a similar relation holds for its free counterpart Uo{t), with Gq{E). 

Note that the analytic structure of G{E) [Gq{E)] in the complex i5-plane is 
symmetric w.r.t. the real E axis, as a reflection of the time-reversal invariance of 
the whole process. Because of the relation |e~*^*| = exp[|i?|t sin(argi5)], we know 
that U{t) [Uo(t)] vanishes for t < if the integration contour C, running from — oo 
to oo along the real E axis, lies a little above all the singularities of G{E) [Go{E)]. 
For the initial- value problem (Cauchy problem), therefore, we have to choose the 
integration contour C in the way shown in Figure 

Let us start our discussion from the case of a system composed of a finite number 
of particles (say, N) put in a finite box (of volume V) . Then H [Hq] has a discrete 
spectrum and, correspondingly, the singularities of G{E) [Go{E)] appear as a series 
of simple poles, running from the lowest-energy point (i.e. the lower bound of the 
spectrum) to infinity on the real E axis. We denote the lowest-energy point in both 
cases by Eg, for simplicity. See Figure |[ 
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Fig. 2. Singularities of G{E) or Go{E) for the case of discrete spectrum and the integration contour 
C for the initial-value problem. 



We can now start the perturbation theory on the basis of the following expan- 
sion 



G{E) = —1— + —1—H'- ^ 



1 ,1 ,1 



E — E — -ff E — H{) 



E — Hq E — Uq E — Hq E — Uq 
^ -H' i H' H' H'- 



The survival probability amplitude (in the Hcisenberg picture) is given by'' 

i 



(26) 



{a\U{t)\a) - — / e~'^'Ga{E)dE, Ga{E) = {a\GiE)\a). (27) 



IC 

We introduce now the plausible assumption that the phases of the off-diagonal 
matrix elements of the interaction Hamiltonian H' are randomized in the infinite 
TV and V Hmit (keeping N/V finite). That is, 

{n\H'\n') have random phases for n ^ n' . (28) 

Note that the discrete energy spectrum becomes continuous in this limit. Conse- 
quently, the above-mentioned simple discrete poles of Gq {E) get closer and closer in 
this limit, and finally merge into a continuous line, distributing over [£"^,00) along 
the real axis, where Eg is the lowest-energy point of the free Hamiltonian. In this 
case, there is a branch point at Eg and a branch cut on [Eg ^00) for G{E) and, in 
particular, for the self-energy part T^a{E) to be introduced shortly. See Figure ||. 

''The quantity Ga{E) is essentially the same as the energy density uJa{E) in (^2[), though the 
meaning of E is different: The latter E stands for the spectrum of the total Hamiltonian H while 
the former E is just an integration variable running from -co to 00 in the Fourier transform, 
irrespectively of the spectrum of H. 
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Fig. 3. Singularities of Ga(_E) for the continuous spectrum and the integration contour for the 
initial-value problem. The branch cut reflects the continuous spectrum of the free Hamiltonian 
and runs along the real E axis. 



Under the above random phase assumption, it is easy to show that Ga{E), which 
is nothing but the Fourier transform of the survival amplitude, is written as 



Ga{E) = 



E-E, 
1 



+ 



1 + 



E-E, 

^a{E) 



{a\{H' + H' 



E-Ea 



+ 



^a{E) 

E-Ea 



E — Hq 

2 



H' + ...)\a) 



1 



E-Ea-^a{E)' 

The "self-energy part" T,a{E) is composed of all the even-order contributions 



s„(E) = s(2)(i;) + sW (£;) + ••• 



where 

^i'HE) 



\{amn)\' 

E — En 



{a\H' i H' i H' ^-^H'\a} 

\{a\H'\n')\^-\{n'\H'\n)\^ 



E-E„. 



E E 



{E - Er,){E - En'Y 



(29) 



(30) 



(31) 



(32) 



and so on. Observe that these quantities correspond to the "proper" self-energy 
parts in field theory (the "improper" parts are subtracted at every step of the 
expansion). In general, we obtain 



^«(^) = E ^ 



E — En 



{E-En'f 



■ (33) 
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In the infinite N and V hmit, of course, the summation in ( |33| ) becomes the following 
integral: 

^a{E) = dE' pa{E')^-^Ka{E') , (34) 

where pa{E') stands for the state number density and Ka{E') for the quantity inside 
the square brackets [• • •] in (|33|). This expression implies that Ga{E) has a branch 
point at Eg and the corresponding branch cut on the real axis, and has no other 
singularities on the first Riemannian sheet, while on the second Riemannian sheet 
possible singularities can appear. See Figure This procedure, based on the above- 
mentioned random phase approximation, is the same as that originally introduced 
in the theory of nuclear reactions It is also very close to the diagonal singularity 
formulated by van Hove 1^1 although we have not yet used his limit concerning the 
time-scale transformation. 

The survival probability amplitude (E9h, with the above Yia{E)^ has been derived 
under the random phase assumption ( |2q ) for many-body systems having a huge but 
finite number of degrees of freedom. However, as is well known, it is an exact relation 
in field-theoretical cases where we have a particle system coupled to fields in free 
space or two (or more) fields coupled to each other. Of course, in such cases the 
number of degrees of freedom is infinite. This is the reason why the quantity Y,a{E) 
has been called the "self-energy part" . Investigation of the analytic properties of the 
amplitude or the relevant 5'-matrix elements on the basis of miantum field theory can 
be traced back to the dampiii? theory initiated by HeitlerO See the review paper 
given by Arnous and Zienau.EJ Therefore, we can discuss the temporal evolution 
of the survival probability amplitude, in the wider framework of field theory and 
many-body problems, on the basis of the expression we have just obtained. 

It is worth pointing out that the familiar Breit-WignerB (or Weisskopf-Wignc 
form does neither match nor follow from the above analysis, since the Breit-Wigner 
spectrum, which is just assumed on a phenomenological basis, extends from —oo 
to oo. Also, notice that if we put E ^ Ea in Tja{E) (which is in general a com- 
plex number), Eq. (|29| ) yields the Weisskopf-Wigncr formula! of the exponential 
decay. Of course, this cannot be justified even in the weak-coupling limit, since the 
exponent should be given by the simple pole of G{E), say E, satisfying 

E^Ea~^a{E), (35) 

which is clearly different from the former, E ^ Ea — T,a{Ea)- 

In order to observe the temporal behavior of the survival probability amplitude 
at longer times, the original contour C, running along the real E axis from — oo to 
oo, is now deformed into a new contour. Recalling again that Ga{E) has a branch 
cut along the real E axis from Eg to oo and a simple pole on the second Riemannian 
sheet, we understand that the {Eg, oo) portion of C is equivalent to the sum of a 
path running just below the real E axis on the second Riemannian sheet from Eg 
to —oo and a circle turning clockwise around the simple pole. The new contour, 
equivalent to the original one, is thus composed of the path C" and the circle on 
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the second sheet, as shown in Figure |^. We therefore obtain from Eqs. (p7|), ( p9| ) 



Fig. 4. The deformed integration contour, composed of a path C' and a circle on the second 
Riemannian sheet. The path C starts from — oo, runs just above the real axis on the first sheet, 
goes down to the second sheet at Eg and runs back to — oo just below the real axis. The circle 
encloses the simple pole clockwise on the second sheet. 

and (H) 

{a\U{t)\a) = ^ I e-''''Ga{E)dE 
J c 

= Ze-*^* + — / e~'^'GaiE)dE, (36) 

where C" is the new contour shown in Figure^ and — [dG~^{E)/dE]^^i^. Since 
the imaginary part of E is negative, the first term, representing the contribution 
from the simple pole on the second Riemannian sheet, yields the exponential decay 
form. As the decay process proceeds further, the exponential term becomes very 
small and the last integral on the contour C", given by 

^ / ' e-'^*[Ga(£; + ie) - Ga{E ~ ie)]dE (37) 

eventually dominates. In the above, the second function in the integrand, Ga{E—ie), 
must be evaluated on the second sheet. Note that [Ga{E + ie) — Ga{E ~ ie)] cx 
\Tja{E + ie) — Tja{E — ie)] cx A'^, where A represents the strength of interaction. The 
power law can be derived from this expression, under appropriate assumptions, 
as an asymptotic formula for very large t. This is in accordance with the general 
theorem (|2^) stated before. The behavior for large t, however, depends on the details 
of Pa{E') and Ka{E'): For this reason, in the next subsection, we shall examine 
the above asymptotic formula in the particular case discussed in the preceding 
subsection. 

It is important here to mention van Hove's limiting procedure concerning the 
time-scale transformation. In spite of their crucial role in the derivation of the closed 
equations (^ and (^ for the survival probability amplitude in the particular case 
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studied in Sec. ||.|l|, the conditions (^) on the interaction Hamiltonian do not seem 
quite natural for systems with many (ideally infinite) degrees of freedom (roughly 
called "macroscopic systems"): The initial state \a) (which is now considered to 
be a direct product of a quantum subsystem and the remaining part of the whole 
system), would be given a rather special status if (^ were to be satisfied. On the 
other hand, the random phase approximation (|2^), introduced at the beginning of 
this subsection, seems more plausible in such a case. PauliEj first derived a master 
equation in the quantum theory of many-body systems along this line of thought. 

Almost three decades later, van Hove showed that the repeated use of the ran- 
dom phase approximation, which is so crucial for the derivation of Pauli's master 
equation, can be avoided by taking the so-called van Hove limilH 

A ^ with X^t kept constant, (38) 

where A represents the strength of the interaction, as mentioned above. He showed 
that the survival amplitude of some suitably chosen initial state decays exponen- 
tially in this limit, provided that the diagonal-singularity assumption holds. This 
assumption, which appears reasonable for systems with many degrees of freedom, 
requires the predominance of the diagonal matrix elements {a\H'H'\a), which turn 
out to be of order N {N representing the number of degrees of freedom of the 
system), over the remaining off-diagonal terms. Observe that in our terminology, 
the limit ( ^Sf ) is equivalent to making the time constant tq, defined in Eq. (|^) of 
Sec. 0, and the integral on C" in (pq) vanish, for the former is ^ 0(A) and the 



latter ~ 0{X^)t {S > 0) [remember Eq. (37) and the subsequent discussion]: In 



conclusion, only the exponential decay, characterized just by the lowest-order con- 
tribution Eq , remains as a purely probabilistic law. By means of this procedure, 
van Hove succeeded in deriving the master equation from the quantum-mechanical 
Schrodinger equation. It is also important to notice that in order to realize the expo- 
nential law, and consequently the master equation, a wave packet must be taken as 
an initial state. It was shownE^ that consideration of wave packets in quantum me- 
chanical scattering processes is equivalent to the extension of scattering amplitudes 
to the complex E plane. The weak-coupling limit (^8|) has also been investigated 
from a mathematical point of view .0^11 

Finally, before closing our general discussion, it is worth showing another method 
to derive temporal behaviors, without resorting to Fourier or Laplace transforms. It 
is not difficult to derive the following differential equation for the survival amplitude 

j^(a\Ujit)\a) 

= - f dh{a\Hi{t)H'j{h)\a){a\Ui{h)\a) 
Jo 

+ [ dh r dt2 [ \t3{a\H',it)H'j{h)H'j{t2)H'iit3)\a)p{a\Ui{t3)\a) 
Jo Jo Jo 

ft pti /•t2 pt^ pt4 

dti / dt2 / dt^ / dt4 / dt5 
Jo Jo Jo Jo 
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x(a|jf;(0i?;(ii)-ff;(t2)i?;(<3)i?;(t4)-ff;(t5)|a)p(a|t//(t5)|a) 



(39) 



on the basis of the iterative solution of the Schrodinger equation 

Ui{t) = l-i [ dt'Hi{t')Ui{t') 
Ja 

^ l-i f dt'H'j{t') + {-if f dh f dt2H'j{ti)H'j{t2) + ■ ■ ■ , (40) 



under the conditions (g) ( "mass renormahzation" ) and (^8|) (random phase approx- 
imation). Observe that due to these conditions, aU terms of odd order give no 
contribution to the survival amplitude. The "proper" matrix elements {a\ ■ ■ ■ \a)p 
in ( |39| ) are defined by (Eat ^ Ea — Eb, etc.) 

{a\H'j{t)H'j{t,)H'j{tM{t,)\a)p 

^ Y,{a\H',{t)H',{t,)\n){nmt2mh)\a) 

= e*''""'^*"*'^e'^"'"^*'"*'^l(a|^f>')n("'l^»l', (41) 

{a\H'j{t)H'jih)H'jit2)H'j{t3)H'jiU)H'jit5)\a)p 

^ E E(«i^/W^/(*i)i")("i^/(^2)i^;(t3)in')(n'ii?;(t4)i/;(i5)ia) 

^ ^ ^ ^ g*-E„„/(t-t5)giS„/„//(ti-t4)giB,.'/„(t2-t3) 

X \{a\H'\n')\^ ■ \{n'\H'\n")\^ ■ \{n"\H'\n)\', (42) 



and so on. It is evident that these terms correspond to the self-energy parts of the 
4th order i:'a\E), 6th order J:''a\E) and so on, in Eq. (^0|), in Fourier space. 

If "memory effects" are neglected so that the amplitudes {a\Ui{ti)\a) appearing 
on the RHS of (^9|) can be evaluated at time t instead of U , the differential equation 
( ^9[ ) is solved to give an exponential form for the survival amplitudeEJ 

{a\Ul{t)\a) ~ e"/o'*^i(*"^i)-''^(^i)+/a''^i/a""'''^^/o""'""' dTs(t-Ti-T2-r,)U{Ti,r2,r3) + - ^ 

(43) 

The functions in the exponent are given by 

/2(ti) = {a\H'e'^^^^-"<'^^^H'\a), (44) 
/4(Ti,r2,r3) = (a|ffV(^"-^«)"iHV(^"-^»)"=i/V(^"-^»)"^i/'|a)p, (45) 

and so on. These expressions allow us to extract some of the typical temporal 
behaviors of the amplitude (except at very large times). To this end, let us recall 
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again that the energy spectrum is asymptoticaUy very dense in many-body systems 
and becomes continuous in field-theoretical cases so that the summations appearing 



in (41 ) and (M), for example, are to be understood as integrations over energy. Here 



we consider the following two casesiEU 

• When the matrix elements of the interaction Hamiltonian H' have a relatively 
wide energy range over which their variations are small, the phase factors 
appearing in the functions /2, /4, ■ ■ • (see ( ^ ) and (^)) oscillate rapidly except 
for very small r's. Therefore, after summing over the intermediate states, the 
functions fi can be considered to be proportional to 6 functions w.r.t. the 
integration variables r's, which results in a linear t behavior 

{a\Ul{t)\a) ~ e'' '^^^M^^'l-Io^ /o°° -^-^ /o°° drsMr,,r,^r,) + -] ^ ^^g^ 

provided the integrals in the exponent converge. Observe that the upper 
limits of integration can be safely taken equal to oo under such conditions 
and therefore the above form becomes a relatively good approximation at 
large t. The time t has been assumed to be much larger than the average 
energy spacing. 

• On the contrary, if the matrix elements have relatively narrow bands, the 
phase factors in the functions /2, /4, ■ ■ • are safely replaced with unity. This 
means that we can replace the arguments of fi in the integrals by zeros and 
arrive at 

(a|;7/(t)|a) ~e-*/^(°)*'+*/*(0'°^o)*'+o(*'). (47) 

Since we have assumed that the time t is small enough so that the above- 
mentioned replacements are allowed, this approximate form is valid at short 
times. We infer that the dominant term at short times is of second order in 
t, in accordance with the argument presented in the preceding section. 

It is worth stressing that the behavior at both short and large times has been 
derived from one and the same expression for the amplitude {a\Uj{t)\a) under dif- 
ferent physical conditions. At the same time, however, we should keep in mind that 
the behavior that really dominates in the particular time region under consideration 
does depend on the relative magnitudes of the Gaussian and the exponential (linear 
t) decays, as well as on the characteristic time constants of the exponents. 

3.3. Reanalyzing the particular case 

In this subsection, we shall reconsider the particular interaction Hamiltonian of 
Sec. and study in detail the analytic properties of the survival amplitude. This 
subsection is mainly written in an educational spirit, and therefore some parts of it 
are a specific version of the general argument presented in Sec. Nevertheless, 
we believe that it is very instructive to explicitly exhibit the temporal behavior at 
long times, and to perform an analysis in terms of Laplace transforms. 
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Let US go back to the definition ( p^ ) of g{s,t) and see how and to which extent 
the exponential law can be realized. Since the amplitude {a\Ui(t)\a) is given by the 
inverse Laplace transform (jll]), we understand that the exponential decay occurs 
only when the complex function g{s, t) has zeros in the left-half complex s plane. In 
order to perform the inverse Laplace transform, the function g has to be analytically 
continued into the left-half complex s plane. By performing the integration over u 
in (|lj) we obtain 

g{s,t) = s + t{a\H'- ^-—^H'\a) 



|2 



Jco -1^0+ 

Here \EQ,r) are eigenstates of iJo, and form a complete orthonormal set, with r being 
quantum numbers describing possible degeneracies. Observe that the integrand has 
a simple pole at Eq = Ea -\- ij. Since t > Q and s is taken to have a positive real 
part (in order to assure the convergence of the integration over t in the Laplace 
transform) , the pole lies above the integration contour Cg which extends along the 
real Eq axis. Therefore in order to extend the function g{s,t) into the left-half 
complex s plane, where the real part of s is negative, the integration contour should 
be deformed in the complex Eq plane so that its relative configuration w.r.t the 
singularity is maintained.^ See Figure ^. This gives us the following change for the 



Fig. 5. Integration contours (a) Co for JR(s) > and (b) Ci for 5R(s) < 0. The former extends 
along the real Eq axis starting from the lowest energy Eg . The pole at Eq = Ea -I- i |- is located 
above the contours Co and Ci . (c) The contour Ci can further be deformed and decomposed into 
the contour Co along the real Eg axis and a circle surrounding the pole. 

last term in (^) 

J Co Ea - Eo+ « J 



JCo Ea- Eq + Ij ^ 



, (49) 
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as the pole moves into the lower-half Eq plane crossing the real Eq axis. Notice 
that the second term of the RHS, which represents the contribution from the simple 
pole now appearing in the lower-half Eq plane, has in general a nonvanishing real 
part. We see that this term just expresses the discontinuity of the LHS over the 
imaginary s axis, thus ensuring the analyticity of the RHS as a whole. 

It is now clear that the analytically continued function l/.g(s, t) has the following 
properties as a function of the complex variable s: 



• Under the assumption that Hq (like H) has a bounded spectrum, a branch cut 
exists along the imaginary s axis, extending to — ioo from the branch point at 
s = it{Ea — Eg), Eg being the lowest value of the spectrum of -ffo- 



• On the first Riemannian sheet, into which the analytic continuation must be 
done without crossing the branch cut, the last term in (|49| ) does not show up. 
The function l/g(s,t) has no singularity and is analytic in this plane. 



• The last term in (49) appears only when the function g(s, t) is continued into 
the second Riemannian sheet through the cut. This term allows the function 
l/g{s,t) to have a simple pole in this plane where 5R(s) < 0. 



These properties of the function l/g{s,t) are enough to understand that the 
survival probability amplitude {a\Ui{t)\a), being nothing but the inverse Laplace 
transform of l/g{s,t), exhibits two distinct behaviors at large times. Observe that 
the integration contour over s in ([ill) can be deformed in the left-half plane as 
in Figure |^(b) . There is a close similarity between the deformed coutours in the 



Fig. 6. Integration contours over s. (a) The original one in (JJj) and (b) the deformed one. The 
latter is composed of a circle surrounding the simple pole at sq on the second Riemannian sheet and 
a path starting from ioo on the second sheet, turning around the branch point at s = it{Ea — Eg) 
and extending back to ioo on the first sheet. 



Fourier (Figure ^) and in the Laplace transforms (Figure H) . Let sq = 



-^t-iSEt 
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(7 > 0) be the zero of g{s, t) in the second Riemannian sheet, so that 
Er\{Eo,r\H'\a)\^ 



-5E 



\Ea-EQ + 5E-q\^ 



dEo 



-2TTY,\{Eo,r\H'\a) 



(50) 



Er\{Eo,r\H'[ 
\Ea-Eo + SE^ 



dEn 



'{Ea-Eo)j:r\{Eo,r\H'\ 



\Ea-Eo + SE- 



-dEn = 0. 



(51) 



Notice that both 7 and SE, which are obtained as solutions of these equations, are 
i-independent constants. Close scrutinyEI shows that in the weak-coupling limit, the 
second term in the square brackets of (|5^) gives a finite contribution and we exactly 
reobtain Fermi's Golden Rule (|l9|) with the right factor. We see again that the 
presence of the second term in the LHS of ( |50| ) is crucial for g{s, t) to have a zero in 
the left-half s plane. This simple-pole contribution to l/g{s^t) provides us with the 
exponential decay form for the survival amplitude, while the integration along the 
imaginary s axis cancels the former to yield the Gaussian behavior at short times 
according to the argument exposed in the first part of Sec. ^. As we shall see, the 
latter contribution is also responsible for the deviation from the exponential law at 
very long times, as expected on the basis of the above-mentioned theorem (p3|). 

To estimate the remaining integration along the imaginary s axis for very large 
times, we chose the integration contour that starts from ioo on the second Rieman- 
nian sheet, goes down along the imaginary s axis, turns around the branch point 
at s = it{Ea ~ Eg), appears into the first Riemannian sheet and runs back to ioo 
again. See Figure |^(b). The contribution along this contour, which we call Xc, is 
explicitly written as^ 



1 



- [{s - itEag) ^ (s - itEag)e-''^'] \ ds. (52) 



27ri J.tEa, V9{s,t) 

Here Eag = Ea — Eg, and the second term is obtained from the first just by replacing 
s — itEag with (s — itEag)e~'^'^^ , thus representing the integration on the second 
Riemannian sheet. With the change of the integration variable s = i{y + tEag), Xc 
is explicitly written as 



Xr 



e "9 

2TTi 









n 




V + tEag 









EJ(i?o,r|i7'| 



H -1 



En — Eq 



-dEn 



(53) 



Let us first calculate the difference between the two integrands. It arises from 
the multi-valuedness of the denominator in the integrand, the relevant part of which 



'^Tho contribution at the branch point is easily shown to vanish. 
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is rewritten, after integration by parts, 

/ p p y — ^-^0 = -^Eq- Eg + ^)y^\{Eo,r\H \ 

+ \n{Eo-Eg + ^)—J2\{Eo,r\H'\a)\'dEo 

= A. (54) 

Notice that the argument of the logarithm is positive because y > 0. Since the 
logarithmic function gains an extra imaginary part — 27ri when its argument turns 
clockwise once around the origin, the above logarithmic function changes into 

\n{Eo -Eg + ^)^ \ii{E„ -Eg + ^)^ 27ri0(f - \Eo - Eg\) (55) 

when y is replaced with ye~^" . Therefore by inserting the above expression for the 
logarithm into formula (|5^), the denominator in the second term of the RHS in ( |53| ) 
is shown to differ from its counterpart in the first term by 

- 27Tit^\{Eg + f ,r|i/'|a)p ee ^27ritB. (56) 

r 

In this derivation, we have assumed that the function J^r I (^o, f P vanishes 
rapidly enough at both boundaries in order to make the integration over Eq con- 
verge: 

x-^ , „ f (£^0 - E„Y for Eo ~ £:„, 6> 0, 

V l(^o)^' forSo^oo, <5'>0. 

We are now in a position to write down the explicit form of Xc- After a few 
manipulations, we arrive at (u ~ y/Eagt) 

^ [Eag{l+u) + Au][Eag{l + u) + Au~2^iBuy ^ ' 

where 

An = / HEn-Eg + Eagu)—Y^\{Eo,r\H'\a)\''dEo, 

J Eg "^0 ^ 

B„ - Y.\^^9+EagU,r\H'\a)\''. (59) 

r 

In the large t limit, i.e., tEag ^ 1, the integration over u can be estimated by the 
behavior of the integrand for small w < mq ~ 1 1 tEag. Recall that the function 
Bu has been assumed to vanish at the boundary E^ = Eg, i.e. at m = 0. As a 
consequence of (|57|), the function Bu behaves, for small u, like 



B^ =Y,\{Eg + EagU,r\H'\a)\' ^EagC^u', J > 0, 

r 



(60) 
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where the first factor adjusts the dimension of B^, while Cq is a positive constant. 
We realize that the function Xc exhibits a power decay at very large times 

Xc ~ 



^ -Ce^^-*^. (61) 

Therefore we have shown that the survival amplitude decays as a function of the 
inverse power at very long times, which eventually takes over the exponential 

behavior. 

To summarize, the survival probability amplitude {a\Ui{t)\a) is shown to have 
the following large-time behavior 

{a\Ui{t)\a) ^ Ze-**-^^^* - Ce^^-*^, (62) 

where Z^^ = {dg{s,t) / ds)\g^_x-iSE- It is very instructive to compare this expres- 
sion with the general formula ( ^ ) of the preceding subsection. It is evident from 
the above derivation that, as was already mentioned, the very large time behavior is 
in fact governed by the state number density pa (E) and the matrix elements Ka (E) . 



See Eq. (|34[ ). (The correspondence B <-> paKa is manifest.) 

The power behavior of quantum survival amplitudes at long times has attracted 
the attention of several investigators in the past. H6hlerE3 showed explicitly a 
power behavior ~ for the Lee model.0 Knight and Milonnill investigated the 
two-level Weisskopf-Wigner model of atomic spontaneous emission under the dipole 
and the rotating-wave approximations. In this model, only two 'essential states' 
are considered, that is, an excited state and no photons present, and the ground 
state and one photon present. They calculated the nondecay amplitude using the 
Laplace transform and found a non-exponential tail ^ at long times. The former 
behavior corresponds to the case S — 1/2 and the latter to 5 = 1 in (|6^). It is worth 
mentioning that these behaviors are in complete agreement with our conclusions, 
since the quantity corresponding to the function i?„ in (p9|) is explicitly given and 



behaves like ^ -^/m (i.e., S = 1/2) in the former case and ^ u {S = 1) in the latter. 



4. Exponential law in a solvable model 



4.1. The AgBr Hamiltonian 

In order to give a concrete idea of how it is possible to obtain the exponential decay 
law for a quantum system in interaction with another quantum system endowed with 
a large number of degrees xif freedom (a "macroscopic system"), we shall consider 
the so-called AgBr model,t3 that is exactly solvable and has played an important 
role in the quantum measurement problem. We present this model in the spirit of 
Havlicek's remark:EJ "Please try to illustrate your assertion on an example which 
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would involve 2x2 matrices only." Our model involves only Pauli matrices and 
some reasonably simple algebra. 

We shall focus our attention on the modified version of the AgBr modcl,Ei that 
is able to take into account energy-exchange processes, and in particular on its 
weak-coupling, macroscopic limit.E3 Our exact calculation will show that the model 
realizes the so-called diagonal singularityEI and can display the occurrence of an 
exponential regime at all times. 

The modified AgBr Hamiltonianll describes the interaction between an ultra- 
relativistic particle Q and a 1-dimensional iV-spin array (D-system). The array 
is a caricature of a linear "photographic emulsion" of AgBr molecules, when one 
identifies the down state of the spin with the undivided molecule and the up state 
with the dissociated molecule (Ag and Br atoms). The particle and each molecule 
interact via a spin-flipping local potential. The total Hamiltonian for the Q + D 
system reads 

H = Ho + H', Ho = Hq + Hd, (1) 

where Hq and Hd, the free Hamiltonians of the Q particle and of the "detector" 
D, respectively, and the interaction Hamiltonian H' are written as^ 



Hq = cp, 



Hn = 



1 ^ 



l + a. 



in) 



H' 



N 



(n) 

cr^ exp 



y^i—x J + a_ exp a; j 



(2) 



where p is the momentum of the Q particle, x its position, V a real potential, 
Xn (n — l,---,Af) the positions of the scatterers in the array {xn > Xn-i) and 
(Tj-'2 the Pauli matrices acting on the nth site._A.n interesting feature of the above 
Hamiltonian, as compared to the original onej£3 is that we are not neglecting the 
energy Hd of the array, namely the energy gap huj between the two states of each 
molecule. This enables us to take into account energy-exchange processes between 
Q and the spin system D. The original Hamiltonian is reobtained in the lo — 
limit. 

The evolution operator in the interaction picture 



Ui{t) 



^iHot/h^~iHt/n 



~i J* H',{t')dt'/h 



(3) 



where H'j[t) is the interaction Hamiltonian in the interaction picture, can be com- 
puted exactly as 



Ui{t) 



N 



T / V{x + ct' 
" Jo 



Xn)dt' 



An) 



exp 



— z— : 
V c 



h.c. 



(4) 



sin this section, we explicitly write the Planck constant h. 
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and a straightforward calculation yields the S'-matrix 

N / vn 

= hm (7/(0 = n • 5(„)=exp -*^rT(").«), (5) 



n=l 



where u — (cos(a;x/c), sin(a;x/c), 0) and Vofl = V{x)dx is assumed to be finite. 
The above expression enables us to compute the "spin-fiip" probability, i.e. the 
probability of dissociating one AgBr molecule: 

If the initial D state is taken to be the ground state |0)jv {N spins down), and 
the initial Q state is a plane wave \p) , the final state reads 



5[^ib,o)A. = ^ .) Hv^r (yT^) ''\pj,j)n, (7) 

3=0 ^ 

where \pj,j)N represents the (spin-symmetrized) state in which Q has energy = 
p — jhoj/c and j spins are up. 

This enables us to compute several interesting quantities, such as the visibility 
of the interference pattern obtained by splitting an incoming Q wave function into 
two branch waves (one of which interacts with D), and the energy "stored" in D 
after the interaction with Q, as well as the fluctuation around the average. The 
final results are 

{Hd)f = qNhuj >nhuj, (8) 



{8Hd)f = ^{{HD-{HD)Ff)F = VmNhuj^V^ncu, 

where f stands for the final state (|7|), p = 1^9, and the trivial trace over the Q 
particle states is suppressed- The arrows show the weak-coupling, macroscopic limit 
N ^ oo, qN = n = finite.E3 All results are exact. It is worth stressing that qN = n 
represents the average number of excited molecules, so that interference, and relative 
energy fluctuations "gradually" disappear as n increases. The limit N oo, qN < 
oo is physically very appealing, in our opinion, because it corresponds to a finite 
energy loss of the Q particle after interacting with the D system. Observe also that 
(^ is a generalized [S'C/(2)] coherent state and becomes a Glauber coherent state 
in the N ^ oo, qN — finite limit. 

4.2. The exponential law 

Our next (and main) task is to study the behavior of the propagator, in order to 
analyze the temporal evolution of the system. Observe that the only nonvanishing 
matrix elements of the interaction Hamiltonian H' are those between the eigenstates 
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of Hq whose spin-quantum numbers differ by one, so that the conditions (Q) in Sec. || 
are not satisfied. It is also important to note that the Q particle state \cp), charac- 
terized by the energy cp, is changed by H' into the state \cpztTiLu), if w ^ 0. We can 
therefore expect a dissipation effect and the appearance of the diagonal singularity, 
which leads to the master equation. Following van Hove's pioneering work,E3 one 
could calculate the propagator perturbatively in this model. However, the solvabil- 
ity of the present model enables us to perform a nonperturbative treatment and 
yields an exact expression for the propagator. Define 

an = an{x, t) = V{x + ct' — Xn)dt' /h, (9) 



which can be viewed as a "tipping angle" of the nth spin if one identifies V with a 
magnetic field and 

4"'(i) = ai"'exp(Tz^S?), (10) 



which satisfy, together with Cg"-' , the SU (2) algebra 



ai"^(J),af (^)J ^ ~ar, (S?), -^"^J = ±2ar(S). (H) 

We can now return to the Schrodinger picture by inverting Eq. (|3|) . We disentangle 
the exponential in Ui by making use of (^fj) and obtain 

N 
n=l 

Notice that the evolution operators (jj) and (|lj), as well as the S'-matrix (H) are 
expressed in a_[actorized form: This is a property of a rather general class of similar 
Hamiltonians.Ej 

We shall now concentrate our attention on the situation in which the Q particle 
is initially located at x' < Xi^ where xi is the position of the first scatterer in the 
linear array, and is moving towards the array with speed c. The spin system is 
initially set in the ground state \Q)n of the free Hamiltonian Hu (all spins down). 
This choice of the ground state is meaningful from a physical point of view, because 
the Q particle is initially outside D. 

The propagator, defined by 

G{x,x',t)= {x\® N{'^\e-'"''^\Q)N®\x'), (13) 



(")^^:^ ^(") 



Ht/h ^~iHot/h ^^—itan{an)<7^_^\3;)^—lncos{an)(T^^^^^—itan{an)(y^_r\^)^ (12) 



is easily calculated from Eq. ( |l2|) and we obtain 

N 

G{x,x',t) = {x\ jv(0|e-'=P*/'' Yl (^e-Mcos[a„(i.t)]}4"'^ \0)n^\x') 

r 

N 

(x\x' + Ct) Yl (e^"^ 



n=l 

N 

{cos[an {x\t)] } 
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N 



= d{x — x' — ct) Y\_ COS an{x' ,t). (14) 

n=l 

Observe that, due to the choice of the free Hamiltonian Hq in (|^), the Q wave 
packet does not disperse, and moves with constant speed c. We place the spin array 
at the far right of the origin {xi > 0) and consider the case where potential V has a 
compact support and the Q particle is initially located at the origin x' — 0, i.e. well 
outside the potential region of D. The above equation shows that the evolution of 
Q occurs only along the path x — ct. Therefore we obtain 

G(x, 0,i) = (5(x- ci) cosa„(i), a„{t) = / V{y ~ Xn)dy/hc. (15) 

This result is exact. Note that 

sin^ a„(cx)) = sin^ ( — ) = q (16) 
\ he J 

is the spin-flip probability (^. We shall now consider the weak-coupling, macro- 
scopic limitEa 



-f-j =0{N-'), (17) 

which is equivalent to the requirement that n — qN be finite. We shall see that the 
limit is equivalent to the van Hove limit, as given by Eq. ( |38| ) of Sec. ^|^. Notice 
that if we set 

Xn=xi + {n-l)A, L ^ XN - xi = {N - 1)A, (18) 

the scaled variable z„ = x^J L can be considered as a continuous one z in the above 
limit, for A/L— s-OasA^— s-cx). Therefore, a summation over n is to be replaced by 
a definite integration 

N 

1 



i2^h{xn)^q-r f h{Lz)dz~n / h{Lz)dz. (19) 

„=1 ^Jxi/L Jxi/L 



This type of integration gives a finite result if the function h is scale invariant, 
because the integration volume is considered to be finite from the physical point of 
view; in fact, the quantities xi/L and xj^/L should be of the order of unity even 



in the L ^ oo limit. It will be shown below [Eq. (20)] that in the present case the 
function h is indeed scale invariant. 

For the sake of simplicity, we shall restrict our attention to the case of J-shaped 
potentials, by setting V{y) = {VQil)6{y). This hypothesis is in fact too restrictive: 
In the following, we shall see that the requirement that V has a compact support 
(local potentials) would suffice. We obtain 

G cx exp f^^lnjcosy" {Von/hc)S{y)dyy\ 



40 Temporal behavior of quantum mechanical systems 



N 



exp [ In {cos [{VQfl/hc)6{ct — Xn)]} 



exp -- 



exp 



x„/L 



d{ct - Lz)dz 



xx/L 
ct — x^ 



(xjv — ct)9{ct — xi) + 9{ct ~ Xn) 



(20) 



where the arrow denotes the weak-coupHng, macroscopic Hmit (p7|). 

This brings about an exponential regime as soon as the interaction starts: In- 
deed, if < ct < Xn, 

to) 



G oc exp 



_c(t 



2L 



(21) 



where to — xi/c is the time at which the Q particle meets the first potential. On the 
other hand, if ct > xn (that corresponds to the case in which Q has gone through 
D and the interaction is over) we have 



(22) 



This could also be obtained directly from ( [l5| ) and is in complete agreement with 
the first formula in (||), because jCp is nothing but the probability that Q goes 
through the spin array and leaves it in the ground state. 

Notice that there is no Gaussian behavior at short times. As was discussed in 
Sec. H and Ref. deviations from the exponential behavior at short times are a 
consequence of the finiteness of the mean energy of the initial state. Observe that 
the mean energy of the position eigenstates is not well-defined in this model. If 
the position eigenstates in Eq. ( p^ ) are substituted with wave packets of size a, the 
mean energy aquires a well-defined meaning. It is shown below that in this case the 
exponential regime is attained a short time after to = xi/c, of the order of a/c. 

Let |0) be a wave packet of size a, initially distributed around the origin x' = 



|0) = / dxC{x)\x) (E) \0)n- 



(23) 

For simplicity, we may choose the following expression for the wave packet 

C{x) = -^eia/2 - IxDe'f"^/'*. (24) 



(One should choose a smooth C function in order to avoid singularities, however, in 
such a case the final expression would be slightly more involved and nothing essential 
would change.) We assume here that the initial wave packet has no overlap with 
the potential region and that the size of the former is much smaller than 

that of the latter 

a/2 < xi, a<^L = XN — Xi. (25) 
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Observe that the mean energy of the initial state is finite, for one obtains 

{0\cp\6) = cpo < oo. (26) 

We calculate the state vectors 



\m) 



|6), 



\Mt)) 



-iHot/h 



|0), 



(27) 



evolved under the action of the total Hamiltonian H and the free Hamiltonian Hq, 
respectively, and define the propagator 



G(t) = {Mt)\m), 



(28) 



which extracts the net effect of the interaction. A straightforward calculation shows 
that the propagator G has the following temporal behavior in the weak-coupling. 



macroscopic limit (17) 



Git) 



1 

.91 W 

92{t) 
93it) 

p-fi/2 



for < ct < xi - a/2, 
for xi — a/2 < ct < xi + a/2, 
for xi + a/2 < ct < xn — a/2, 
for xn — a/2 < ct < xn + a/2, 
for Xn + a/2 < ct, 



(29) 



where the functions gi(t) (i — I, 2, 3) are given by 



9iit) = 
92{t) = 

53 (i) = 



1 + ^ 
na 

2Lr 

na - 



— exp 

-na/2L 



-a/2) 



1, 

[ct — Xl 

a 



a/2). 



exp 



/ n , 



Xl 



2L r 
1 + — 

na 



exp 



ict- 

2L^ 



Xn 



-a/2) 
a/2) 



1 



-{ct 
a 



Xn ■ 



(30) 
a/2) 



See Figure 7. 

We realize that only the function 52 (i) exhibits the exponential form and that in 
the time intervals x\/c—a/2c < t < xi/c+a/2c and xn / c~a/2c < t < xn /c+a/2c, 
the propagator G does not follow an exponential behavior. In the first interval the 
wave packet has not yet completely entered the potential region and the interaction 
between the Q particle and the spin system D has not yet attained its full strength. 
Analogously, in the second region, the wave packet is going out of the array. Notice 
that in this model the exponential law holds if the whole wave packet is merged in 
the potential region so that the interaction keeps its full strength: See 52 (i) in (|30|). 
This suggests an equivalent, maybe more "suggestive" interpretation: one could say 
that in the first region xi/c — a /2c < t < xi/ c + a /2c, the "initial state" has not 
been "prepared", yet. (By "initial", in the above, we do not mean an eigenstate 
of the free Hamiltonian, but rather, intuitively, the state that will undergo the 
exponential decay.) This alternative viewpoint sheds some light on the importance 
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Fig. 7. Temporal behavior of the propagator G{t) defined in See 



of what is usually referred to as "state preparation" in quantum mechanics: It is first 
necessary to prepare the right initial state, in order to observe a certain temporal 
(e.g. exponential) development of this state. 

The region xi/c~ a/2c < t < xi/c + a/2c may be viewed as a possible residuum 
of the short-time Gaussian behavior. In fact, in this region, the function gi{t) 
behaves like 

^(rf-xi+a/2)2, (31) 

and a similar behavior of gsit) is found in the final region of the linear array. It 
is worth stressing that the temporal behaviors ohtaiued in (|l|) and in ( p^ ) are in 
complete agreement with some general theorems .ElO'ESl 

The absence of the power behavior at very long times in this model may be traced 
back to the special character of our free Hamiltonian Hq in (||): The spectrum of 
the operator cp is not bounded from below, so that the total Hamiltonian H does 
not possess a lower-bounded spectrum. Of course we can safely consider, on physical 
grounds, only those incident Q particles that have enough (positive) energy to go 
through the spin system, in order to assure the positivity of the energy of the total 
system. Nevertheless, it should be recalled that the boundedness of spectrum of 
the total Hamiltonian is the only condition required for the Paley- Wiener theorem 
( p3| ) in Sec. |^. Our model system is out of the range of applicability of the theorem 
and therefore its exponential behavior at very long times is not in contradiction 
with it. Similarly, no ground state of the free Hamiltonian Hq in (|^) exists owing 
to the presence of the cp operator and the analysis developed in Sec. ^.|| is not 
applicable to this case either. In this respect, it would be interesting to consider 
the lower-bounded Hamiltonian Hq — c\p\, but, unfortunately, this operator raises 
serious problems that are closelvrelated to the early atttempts at defining a time 
operator in quantum mechanics .O 

One might suspect that the approximation of (5-shaped potential is playing an 
important role in the derivation of the exponential law. As a matter of fact, this is 
not the case: A detailed calculation, making use of square potentials of strenght Vq 
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and width fi, centered at each a;„, yields, for xi + ^ < ct < xn ~ ^, 

Gcxexp(^-^^ + — j. (32) 

In this case, the exponential regime is attained a short time after to, of the order of 
the width of the potential V, which, in the present model, can be made arbitrarily 
small. Also in the present context, the region t ~ tg + 0{^l/c) may be viewed 
as a residuum of the short-time Gaussian behavior. However, it is clear from Eq. 
( p2[ ) that once the particle has completely penetrated into the potential region, the 
exponential law holds exactly in the weak-coupling, macroscopic limit. This is, in 
our opinion, one of the most interesting features of the AgBr Hamiltonian. 

It is very interesting to bring to light the profound link between the weak- 
coupling, macroscopic limit qN ^ n — finite considered in this model and van 
Hove's "A^T" limit.E3 First, it is important to note that in the weak-coupling, 
macroscopic limit, van Hove's diagonal singularity naturally appears in the present 
model. It is easy to check that for each diagonal matrix element of H'^ , there are 
N intermediate-state contributions: Indeed, for example 

N 

(0, . . . , 0|ff'2|0, . . . , 0) = ^ 1(0, . . . , 0|il'|0, . . . , 0, !(,), 0, . . . , 0) p. (33) 

On the other hand, at most two states can contribute to each off-diagonal matrix 
element of H'^ . This ensures that only the diagonal matrix elements are kept in the 
weak-coupling, macroscopic limit, N ~f oo with qN < oo, which is the realization 
of diagonal singularity in our model. 

The link between the weak-coupling, macroscopic limit qN = n = finite con- 
sidered above and van Hove's "A^T" limit can be also intuitively justified as fol- 
lows: The free part of the particle Hamiltonian is Hq = cp, so that the par- 
ticle travels with constant speed c, and interacts with the detector during the 
total time T = L/c ~ NA/c. Since the coupling constant A oc VqH, one gets 
A^T cx (Vbil)^A^A/c cx qN. Notice that the "lattice spacing" A, the inverse of 
which corresponds to a density in our 1-dimensional model, can be kept finite in 
the limit. (In such a case, we have to express everything in terms of scaled variables, 
that is, T = t/L, zi and zjy, introduced just after Eq. (|l^ ) and C = a/L, where a is 
the size of the wave packet.) 

The role played by the energy gap cu need also be clarified. First of all, u 
plays a very important role to guarantee the consistency of the physical framework: 
If = 0, all spin states would be energetically degenerated and the choice of 
the (T3-diagonal representation would be quite arbitrary. In other words, only a 
nonvanishing lo (or Hjj) logically enables us to use the eigenstates of in order 
to evaluate the relevant matrix elements. Although lu does not appear in our final 
results (^ ) or (^9|), it certainly does in other propagators, involving an initial state 
of the spin system different from |0)7v- 
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Finally, we would like to shortly comment on other possible causes for the oc- 
currence of the exponential behavior displayed by our model. This is a delicate 



problem. Our analysis suggests that the exponential behavior (21) is mainly due to 
the locality of the potentials V and the factorized form of the evolution operator 
(|l2|), which shows that the interactions between Q and adjacent spins of the array 
are independent, and the evolution "starts anew" at every step. This suggest the 
occurrence of a sort of Markovian process, which would be in agreement with the 
purely dissipative behavior . A preliminary calculation shows that the operator 

W{t) = cm - {cp{t)). (34) 

where p{t) is in the Heisenberg picture and the average is taken over a wave packet 
of Q and the ground state of Z?, has several interesting features: In particular, some 
properties of W{t) closely resemble those of a Wiener process. In this context, the 
connection between the exponential "probability dissipation" (|2^) and the (practi- 
cally irreversible) energy-exchange between the Q particle and the "environmentii 
-D is a very open problem and should be investigated in detail. Leggett's remark,E3 
about the central relevance of the problem of dissipation to decoherence and quan- 
tum measurement theory make the above topic very interesting: There is, in our 
opinion, a profound link between decoherence, irreversibility, "probability dissipa- 
tion" and genuine energy dissipation. This link is not yet fully understood. Work 
is in progress to clarify the above issues. 



5. Conclusions and additional comments 

We have analyzed the temporal behavior of quantum mechanical systems at all 
times, by taking into account seminal as well as more recent results. As we have 
seen, the temporal development of the survival probability of a suitable initial state 
undergoes an evolution that can be roughly decomposed into three regions: A 
Gaussian region at short times, an approximately exponential region at intermediate 
times and a power region at very long times. 

The short-time Gaussian behavior is ascribable to a very general property of 
the temporal development engendered by the Schrodinger equation. As we have 
seen, this leads to the quantum Zeno effect and has very interesting spin-offs on 
the quantum measurement problem. In particular, we have argued that the QZE 
is by no means a proof in support of the Copenhagen interpretation, based on von 
Neumann's projecion postulate: Indeed, a purely dynamical explanation of the QZE 
was put forward, completely based on unitary evolutions. We have also discussed 
the important difference between the QZE and the QZP, and have argued that the 
latter is in contradiction with Heisenberg's uncertainty principle. 

One of our primary concerns has been to understand under which physical con- 
ditions we can observe an exponential behavior, since, as was already stressed, no 
deviation from the exponential law has ever been observedJlj The exponential form 
is representative of a totally incoherent nature of the system under consideration. 
Remember that the notion of coherence is one of the essential features of the quan- 
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turn world, where the evolution (of the total system) is always unitary. When we 
focus our attention on one part (subsystem) of the total system, if the latter is 
endowed with a huge number of degrees of freedom and can be treated as a kind of 
reservoir, we can expect the former to behave in a stochastic manner. The interac- 
tion between the subsystem and the remaining part of the total system may, under 
such circumstances, cause dissipation on the former. It is then natural to expect 
the appearance of irreversible phenomena, governed by master equations, for such 
subsystems. This is, however, by no means a trivial matter: Because quantum dy- 
namics, i.e. the Schrodinger equation, is symmetric in time, we obviously expect the 
survival probability to go back to unity after a certain elapse of time, which is of 
the order of Poincare time. However, a purely exponential decay of the amplitude 
would actually imply an actual ( "eternal" ) loss of probability, that would never be 
recovered. Rigorously speaking, this is at variance with the underlying unitarity of 
the temporal evolution. An appropriate limiting procedure is required in order to 
explicitly derive the exponential law, as was repeatedly remarked. 

Since it is mathematically true that quantum systems never behave exponentially 
at very long times, it is clear that in order to realize the exponential law at long 
times on the basis of quantum theory, some additional manipulation is necessary. 
By considering the exponential law as a manifestation of the classical character of 
the system, we can recognize the importance of such a notion as the transition from 
a quantum to a classical regime. 

Fermi's Golden RuleJ^ according to which we obtain the decay constant ( p^ in 
Sec. ^ is a lowest-order approximation. We believe that the quantum mechanical 
lowest-order perturbation theory cannot be considered completely satisfactory, in 
this context, unless it is appropriately supported by some additional limiting pro- 
cedure that accounts for the macroscopic nature of the total system with which the 
decaying system is interacting. 

We have also seen that van Hove's procedure, based on the weak-coupling, 
macroscopic limit considered in Sec. ^, yields a satisfactory answer to many of 
the above problems. In particular, his procedure dispenses with Pauli's random 
phase assumptions and yields an exponential behavior at all times, by cancelling 
both the short-time Gaussian region and the long-time power tail. Of course, one 
is led to wonder about the meaning of the "A^T" limit, and in particular about the 
very reasons why Nature should follow this prescription. Even though no clear-cut 
answer can be given to this very deep question, we believe that there is much more 
to be understood, on this issue. 

What we are actually looking for is a consistent way of deriving the classical 
nature (loss of coherence) from quantum theory in a certain limit which should 
properly characterize the classical and macroscopic nature. At this poin t one may 
see a close connection with the quantum measurement problem,E2I^E3c3 where the 
realization of the classical nature of experimental results, namely the appearance of 
exclusive events after a measurement, has been one of the central issues of contro- 
versy since the birth of quantum mechanics. 
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As a general rule, explicit examples can be often more useful than general ab- 
stract considerations. As Berry recently put it. "Only wimps specialize in the 
general case. Real scientists pursue examples" In the attempt to follow the same, 
general philosophy, we have endeavoured to clarify the above issues by considering 
a particular solvable model, known as modifed AgBr Hamiltonian: This model en- 
joys several interesting features: It is solvable and relatively simple. It also yields 
very useful insights into general topics like quantum measurements and dissipation. 
Finally, it discloses the occurrence of van Hove's limit in a very direct way. 

As we have seen, if the initial state of the spin array in the AgBr model is chosen 
to be the ground state of the free Hamiltonian, in which all spins are down, the 
survival probability of such a state, once the interaction with an incoming external 
particle has been switched on, behaves in a purely exponential way. The causes 
of this behaviour, as well as the physical approximations they reflect, have been 
analyzed in detail, also in the light of general theorems. 

In conclusion, we believe that further understanding in all the above-mentioned 
topics can be achieved by analyzing the general structure of quantum mechanical 
evolutions and, at the same time, by looking at interesting models and particular 
limiting procedures. 
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